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Abstract 



We show that, in any conformal field theory, the weights of all bulk primary 
fields that couple to N 02 l fields on the boundary are given by the spectrum 



of an iV-particle Calogero-Sutherland model. The corresponding correlation 
function is simply related to the iV-particle wave function. Applications are 



CO 

discussed to the minimal models and the non- unitary O(re) model. 



The quantum Calogero-Sutherland (C-S) model has proved to be ubiquitous in theoretical 
physics. It has arisen in various ways in conformal field theory (CFT) in the past [1]. In this 
note, we point out a very direct connection. This was originally discovered [2] in the course of 
developing a multiparticle generalisation of Schramm-Loewner Evolution (SLE) [3], which, 
largely through the work of Lawler, Schramm and Werner (LSW) [4] , has recently enlarged 
our perspective on conformally invariant random processes. However, the connection to the 
C-S model may be derived independently from SLE, using the basic principles of CFT, and 
it will now be presented in this way. 

The set-up is as follows: suppose we have a CFT in the interior of the unit disc \z\ < 1, 
with a conformal boundary condition, and consider in particular the correlation function 

(0(e^)...0(e^)$(O)) = ^ 1 ,...,^|$) (1) 
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of TV boundary fields with a single primary bulk field $ at the origin. (Of course, any 
correlation function in a simply connected region with N boundary fields and a single bulk 
field at an interior point may be related to this by a conformal mapping.) In the second ex- 
pression we have written this correlation function in the operator formulation of CFT, using 
radial quantisation: here |<&) is a highest weight state of the holomorphic and antiholomor- 
phic Virasoro algebras, and |6*i, . . . , 6^) is a state given by the action of boundary operators 
4>(e z9j ) on the boundary state corresponding to the given conformal boundary condition. 
These states lie in an iV-dimensional subspace of the full Hilbert space of the CFT. 

Let us now suppose that is a primary field which is degenerate at level 2: it is a 02, i 
(or 12 ) field in the Kac classification. As shown many years ago by Belavin, Polyakov and 
Zamolodchikov [5], this implies that correlation functions such as (1) satisfy second-order 
differential equations. In this case we shall show that these imply the C-S equation. 

First fix some notation: parametrise the central charge by c = 1 — 6(4 — k) 2 /Ak, so that 
the boundary scaling dimension of is h 2> i = (6 — k)/2k, and the null vector condition is 
(L_ 2 — (k/4)L 2 1 )|02,i) = 0. Define the iV-particle C-S hamiltonian with parameter (5 by 

h im - 1 f 02 + 13(13 ~ 2) v 1 m 

and the free fermion wave function 

^(01,...,^)= I] (e*'-e**) (3) 

l<j<k<N 

Then (subject to suitable boundary conditions, see later) the ground state wave function of 
H N {(3) is l^l^ 2 with energy ({3/2) 2 E% where E% = ±N(N 2 - 1). 

Our main result is that if the correlation function (1) is non-vanishing, then the highest 
weight (bulk scaling dimension) of $ is given by 

x* = h» + h* = ^An(8/k) - ^-E* + + ^ (4) 

where Ajv(/3) is some eigenvalue of H^{f3). Moreover, the correlation function (1) is (up to 
a normalisation) equal to the corresponding eigenfunction, divided by |^/at| 2 / k . 



Consider the infinitesimal conformal transformation z 



z + a(z), where a(z) 



T,j=ibj(Xj(z), with 




z + 



(5) 



= —z 



z — e 3 



where the bj are infinitesimal parameters, initially chosen to be arbitrary. Note that this 
preserves the unit circle with the points {e l6j } removed, but near the origin it acts as a 
pure dilatation. It may be implemented by inserting into the correlation function (1) a 
factor J c T(z)a(z)dz/2ni — f c T(z)a(z)dz/2iri, where T and T are the holomorphic and 
antiholomorphic components of the stress tensor, and C is any contour in \z\ < 1 which 
encircles the origin once in a counter-clockwise sense. 

The effect of this insertion may be evaluated in two ways: first by shrinking the contour 
towards the origin and observing that, as z — > 0, a(z) ~ z + 0(z 2 ), and then using the 
operator product expansion (OPE) T(z)&(0) = h$&(0)/z 2 + 0(z^ 1 ), together with a similar 
antiholomorphic expression. This recovers the original correlation function multiplied by 

+ h$) Y,j bj- It is important for this that $ is primary, so that the higher order terms in 
the expansion of a(z) do not contribute. 

The other way is to distort C so as to lie along the unit circle, with small semicir- 
cles excluding the points {e j }. The contributions from the parts of the contour on the 
circle vanish by virtue of the conformal boundary condition [6]. That from the semi- 
circle around e ld] may be evaluated by writing z = e* ej+l< >, where ( is a local coordi- 
nate whose imaginary part is zero along the boundary. Expanding aij in powers of (, 
we find, after a little algebra, aj = bj(2/( — £/6 + 0(( 2 )). Using the OPE with the 
stress tensor again, the effect of the infinitesimal transformation ctj on (j)(e l6j ) is to gen- 
erate — 6j(2L_ 2 — \Lo)<j)(e ldi ) = —bj((K/2)(d/d6j) 2 — \h2^)4>{e tdj ) (the minus sign is because 
C wraps around e l9j clockwise.) On the other hand, aj is regular at the points e l6k with 
k ^ j, so that 0(e iefc ) -> (1 + 6 j ^(e iefc )) /l2 - 1 0(e^ + b jaj (e i9k )). 

Putting together these contributions, the effect of otj on the boundary state is equivalent 

to 
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k d 2 1 _ / e k - e, a . e k - e 3 1 ' 

■ — 7^7 H ^2 1 — > c °t T^T" + 1 COt -flo 1 177- . . /^o 1 

2^] 6 2,1 ^. V 2 <9£ fc 2 2>1 2sin 2 (0 fe -0 j )/2 2 '\ 



(6) 



Now sum over j: the penultimate terms in the above sum to something proportional to 
J2j Hk^jipj + b k ) cot(9 k — 0j)/2, which vanishes if we now take bj = b k for each pair (j, k). 
The generator of the transformation, acting in the subspace of boundary states, is therefore 

G = "2 S «J - Sgcot — ^ + (Eg 28in% _ W2 T J fc « P> 
The first two terms can be recognised as a similarity transform of H n (4/k), up to a constant: 
|* K p G |*„|- 2 /* = K [ffK(4/K) _ (2/K)X ] + ^ g + f ) fe,, (8) 

The penultimate term then combines with the potential term in Hn(A/k) to give a Calogero- 
Sutherland hamiltonian at a shifted value of /3, equal to 8//c: 

IWjvI^GI^I- 2 /" = - (4/«)£* + ^/i 2>1 (9) 

The eigenvalues and eigenvectors of G and H n (8/k) are thus simply related, and since the 
left eigenvalues of G are N(h^, + h^), we get (4), except for the last term. This arises because 
the same insertion must be made in the partition function, which transforms non-trivially 
owing to the presence of a non-zero trace (6) of the stress tensor at the curved boundary 
[7]. However, a clearer derivation of this term may be found by considering the conformally 
equivalent geometry of a semi-infinite cylinder parametrised by the complex variable \nz. 
The boundary is now no longer curved, but the dilatation operator becomes the generator 
of translations along the cylinder, which is [8] Lq + Lq — c/12. 

We now discuss some examples and applications of the general result. The C-S hamil- 
tonian Hjsr commutes with the total momentum P = —iJ2j(d/d9j), whose eigenvalues cor- 
respond to the spin h$ — h$ of the bulk primary field. The eigenvalue equation admits 
two possible boundary conditions on the wave function ip as a given pair (j, k) of particles 
approach each other: ip oc \9j — ^| 7 , with 7 = (3/2 ('fermionic') or 1 — (3/2 ('bosonic'). 
These correspond to the values allowed by the BPZ fusion rules [5]: the correlation function 
behaves as \9j — 9 k \ 1 ~ 2 l K ) which is consistent with the OPE 



<hM) ■ <t>2M) ~ & - ^r 2 " 2 ' 1 1 + & - k \ h3 - 1 - 2h2 ' 1 03,1 (io) 

with /i 2 ,i = (6 — k)/2k and ft^i = (8 — k)/k. 

The eigenvalues of H^{f3) in the fermionic case are well-known [9]: they have the form 
A = \ J2f=i kj, where the allowed values of the quasiparticle momenta kj satisfy Y,j kj = P 
and kj + i — kj = \fi + Pj, with pj a non-negative integer. The ground state in this sector has 
all the pj = 0, and eigenvalue ((3/2) 2 Efj, which, after a little algebra, leads to a weight 

^-^-ii^l! (11) 
Xn ~ 2k 8k (Uj 

of the corresponding bulk field. 



The 0{n) model. 

The most immediate example of such a bulk field is in the non-unitary, non-minimal 
CFT that is supposed to represent the scaling limit of the O(n) model [10] with n G [—2, 2]. 
This may be realised as a gas of non-intersecting closed loops and open curves [11], in which 
open curves ending on the boundary are known to be described by 2 ,i fields [6]. The above 
result for x { N then agrees with the known value [12] for the bulk TV-leg field. That is, the 
correlation function (1) is proportional to the probability that N non-intersecting curves 
connect the origin to the points e l6j on the boundary. 

For N — 1 and P ^ we find a weight = x[ + (k/2)P 2 with P an integer. The 
physical interpretation of these new primary fields is in terms of winding number states: 
each curve linking the origin and the boundary with winding number \ is weighted by e tFx . 

For iV = 2, in addition to the winding states with P^O, there are new spinless primary 
fields corresponding to k 2 = —k\ = (3/ A + p, with weights x$ = x\ + |/tp 2 + 2p, where 
p is a positive integer: these correspond in the O(n) model to excited modes of the pair 
of curves, confined by their mutual repulsion. In the case of a finite cylinder, or annulus, 
their weights give the exponents of correction terms in the correlation function. There are 
further N = 2 primary fields corresponding to purely bosonic boundary conditions, with 
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7 = 1- (3/2. These correspond to A 2 = (1 - [3/2) 2 E%, which gives x\ = for all k. This 
is consistent with the OPE (10): the two boundary fields are fusing to the identity on the 
boundary, which then couples to the identity field in the bulk. However, once again there 
are excited states in this sector, which correspond to possible new bulk primary fields with 
weights x$ = ^np(p + 1) — 2p. These two types of boundary condition may be understood 
within the O(n) model as follows: each boundary field 2 ,i and its attached curve carry an 
O(n) vector index. If the two labels are different, the curves cannot join, and the fusion is 
into a 03 ; i field transforming according to a tensor representation of O(n). However, if the 
labels are the same, the curves can join up before reaching the origin. The fusion in this 
case is into the identity field, and the fact that the leading coupling is now to a bulk primary 
field with x$> = means that the probability of the curves joining, and not therefore passing 
thorugh the origin, is unity. The weights corresponding to the excited states, with p > 1, 
then give the exponents of correction terms to this for an annulus. p — 1 corresponds to the 
bulk energy density field of the O(n) model. 

However, still for N = 2, there are other possible 'mixed' boundary conditions which 
are fermionic as 9 2 — Q\ — > 0+, and bosonic as 9 2 — Q\ — > 2-k— . The ground state in this 
sector has energy A 2 = jq, which gives x$ = (3k — 8) (8 — re)/32/c. This is the exponent 
determining the relative probability that a curve, whose ends are attached at nearby points 
on the boundary, should enclose the origin or not. For k — 6, this is the 'one-arm' exponent 
of percolation (related to the probability that the origin lies in a cluster which touches the 
boundary), as computed by LSW [13]. Once again, there are excited states in this sector 
whose energies give corrections to scaling. 



Minimal models. 

These correspond to rational k = 4k/ k', and the allowed values of the weights of scalar 
bulk primary fields are given by the Kac formula x r>s = ((4r — ks) 2 — (4 — k) 2 )/8k with 
1 < r < k — 1 and 1 < s < k' — 1. For these to agree with (4) imposes a severe constraint. If 
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it cannot be satisfied, it implies that the correlation function (1) must vanish. For example, 
the weight of the 1-leg field corresponds formally to x 1 / 2j o, but for this to appear in the table 
of allowed values it is necessary that k is odd, with r = (k — 1)/2, and k! even, with s = k'/2. 
This is of course consistent with the fusion rules of boundary conformal field theory. For 
N = 2, with fermionic boundary conditions, no solution for (r, s) in the allowed range is 
possible, indicating that, in a minimal model, two 02, 1 fields on the boundary can couple to 
the bulk only through fusion into the identity. In that case, coupling is allowed, as long as 
p < [k/2] — 1, because the corresponding allowed weight is x$ = xi,2 P +i- 

Comparison with multiple SLEs. 

Although the above CFT arguments are self-contained, it is instructive to compare them 
with those of Ref. [2]. There, a multi-particle generalisation of SLE was proposed in which 
the N curves connecting the boundary with the origin are 'grown' dynamically, starting from 
the boundary at time t — and reaching the origin as t — > oo. This process is described 
in terms of the evolution of the conformal mapping gt(z) which sends the simply connected 
region not yet excluded by the curves into the whole unit disc. This turns out to satisfy 
dg t = J2jbj®j(gt), with otj having the same form as in (5). In this picture, however, the 9j 
become functions of t, evolving, if we take bj = dt for all j, according to Dyson's brownian 
motion [14]: 

dBj = J2 cot((0j - 9 k )/2) dt + dBj(t) (12) 

Mi 

The terms on the right hand side correspond to a mutual repulsion and a stochastic noise. 
For N — 1 it is known [3] that this process (radial SLE) gives the correct measure on the 
continuum limit of a single curve. The corresponding Fokker-Planck equation for the joint 
probability distribution P({9j};t) has the form dP/dt = CP, and it was argued in Ref. [2] 
that the asymptotic equilibrium solution of this equation, satisfying CP eq = 0, should give 
the distribution of the points on the boundary in an equilibrium 2d critical system such as 
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the 0(n) model. This gives the result |\1/7v| 4//k , in contradiction with the CFT prediction 

|^v| 2/K - 

This discrepancy may be traced to the form of 

The first term comes from averaging over the white noise dBj, and the second from the 
repulsion. Comparing this with the expression (7) for the conformal generator G, we see 
that the last term, proportional to h 2 ,i, is absent. This arises, in the CFT calculation, from 
the transformation of (f>(e l9k ) under the local scale transformation induced by the conformal 
mapping otj, but it is missing in the multiple SLE approach. Indeed, there is no obvious 
way of modifying (12) so as to incorporate such a term. 

However, the argument in Ref. [2] that the joint distribution of the boundary points {9j} 
in a critical system should be given by the equilibrium distribution of the process (12) was 
also based on the assumption that the measure on the N curves was strictly conformally 
invariant under g t . For the case of a single curve with each end on the boundary of a simple 
connected region, this is known not to be the case in general - rather, it is conformally 
covariant. [4] The covariance factor is just the product of the local scale transformations at 
each end, raised to the power h 2 ,i- When the arguments of Ref. [2] are modified to take 
into account these factors, the result agrees precisely with that of the CFT argument given 
earlier. The fact that SLE K is now associated with H N (f3) with (3 = 8/k is much more 
satisfactory. For example, it suggests that the duality of SLE [15] under k — > 16/ k, is related 
to the known duality of the C-S hamiltonian [16] under (5 — > 4//3. 

To summarise, we have shown that the quantum Calogero-Sutherland model arises in 
a very simple way in bulk-boundary conformal field theory. The full spectrum is realised 
in the non-unitary CFT of the O(n) model, and it predicts the scaling dimensions of new 
primary fields in that theory. In minimal models, it places severe restrictions on which bulk 
fields can couple to the boundary. In general, because of the Galilean invariance of the C-S 
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hamiltonian, any scalar primary bulk operator is associated with a tower of other primaries 
of spin s, with the differences in conformal weights proportional to s 2 . Thus, although the 
spectrum of Virasoro descendants in a CFT is relativistic, the spectrum of these primaries 
has a non-relativistic form. 

Although we have considered only the case of 2 ,i fields on the boundary, our arguments 
can be generalised to include other boundary fields corresponding to degenerate Virasoro 
representations. These will lead to higher-order differential operators. Similar generalisa- 
tions to WZWN models are also possible. 

Recently Bauer and Bernard have extended their analysis of CFT as a probe of SLE [17] 
to the radial case [18]. Some of their results overlap with ours. 

The author acknowledges both helpful and critical conversations with D. Bernard, 
G. Lawler and W. Werner. This work was supported in part by the EPSRC under Grant 
GR/J78327. 
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